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Introduction. 

In a paper published in ' Acta Mathematiea/ vol. 16, M. Folie announces the fact 
that the latitude of places on the earth's surface is undergoing periodic changes in a 
period considerably in excess of that which theory has hitherto been supposed to 
require. This result has been confirmed in a remarkable manner by Dr. S. C. 
Chandler, in America (vide ' Astronomical Journal/ vols. 11, 12), who, as the result of 
an exhaustive examination of almost all the available records of latitude observations 
for the last half-century, has assigned 427 days as the true period in which the 
changes are taking place. 

The old theory, based on the assumption that the earth was rigid throughout, led 
to a period of 305 days, and M. Folie proposes to account for the extension of this 
period by attributing a certain amount of freedom to the internal portions of the 
earth. The earth he supposes to be composed of " a solid shell moving more or less 
freely on a nucleus consisting of fluid at least at its surface. " The argument advanced 
by M. Folie in favour of this constitution of the earth, namely, the independence of 
the motions of the shell and the nucleus, appeared to me to be unsatisfactory, and I 
therefore proposed to myself to test the validity of it by examining a particular case 
which lent itself to mathematical analysis, namely, that in which the internal surface 
of the shell is ellipsoidal and the nucleus consists entirely of homogeneous fluid. 

The principal axes of the shell and of the cavity occupied by fluid are assumed to 
be coincident, and the oscillations are considered about a state of steady motion in 
which the axis of rotation coincides with one of these axes. It is clear that a steady 
motion will be possible in this case, and that such a motion will be secularly stable 
in the event of the axis of rotation being the axis of greatest moment for both the 
shell and the cavity. 

The problem was originally treated by the analysis used by Poincake in his 
memoir on the stability of the fluid ellipsoid with a free surface (' Acta Mathematiea/ 

25.7.95. 
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vol. 7). This analysis reduces the determination of the motion of the fluid to the 
problem of finding a single function i//, subject to certain boundary conditions, which 
in our case take a very simple form. In the case where the surface of the fluid is 
ellipsoidal, it is found that, when the system, is oscillating in one of its normal modes, 
\jj will be expressible as the sum of a series of Lame products of a single order n only. 

When n is different from 2, the types of oscillation are such that no disturbance of 
the shell is involved, and a period equation for the oscillations of the fluid may be 
deduced in a manner similar to that riven by Poincare. 

The types of oscillation corresponding to n = 2 demand exceptional treatment, in 
consequence of the motion communicated to the shell when they exist. The fluid 
motion, however, is found to be such that the molecular rotation is everywhere the 
same. Mr. Bryan has suggested to me that this circumstance may be made use of 
in order to treat the oscillations which involve motions of the shell by a simple 
analysis previously employed by Greenhill ('Proc. Camb. Phil. Soc./ vol. 4, p. 4) 
which does not involve Lame functions. To facilitate the reading of the paper, the 
results are first deduced by this method, and the Lame analysis by which they were 
originally obtained is reserved for an appendix. 

The oscillations under consideration are found to be of two types. One of these 
corresponds to an oscillation previously discussed by Hopkins in his " Researches in 
Physical Geology " ('Phil. Trans./ .1839). This exists only in consequence of the 
contained fluid, and in it the oscillations of the shell are similar in character to the 
" forced " nutations of the earth produced by the action of the sua and moon. In 
the other type the motion of the shell is closely analogous to the motion of a rigid 
body when slightly disturbed from a motion of pure rotation about a principal axis, 
and, in fact, identifies itself with such an oscillation in the event of the inertia of the 
fluid becoming negligible. 

On applying the problem to the case of the earth, the latter mode is that on which 
the variations of latitude depend. The period^ however, is found to be shorter than 
it would be if the fluid were solidified, and thus, in this particular case, M. Folie's 
results are contradicted. It appears to me to be highly probable that any such 
freedom in the interior of the earth as that supposed by M. Folie, provided the 
surface does not undergo deformation, would have the effect of reducing, instead of 
extending, the period, and the true explanation of the phenomenon is probably that 
given by Newcomb (' Monthly Notices of the Royal Astronomical Society,' March, 
1892), who shows that the elasticity of the earth, as a whole, would have the effect 
of prolonging the period. 

§ 1 . The Period Equation. 

Let us refer to rectangular axes coincident with the principal axes of the ellipsoidal 
ca vi ty. 
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Let a, ft, y be the principal semi-axes of this cavity ; A, B, C the principal 
moments of inertia of the shell. 

Suppose the motion of the fluid at any instant consists of a rigid body rotation 
with angular velocities f, 77, £ about the axes of the ellipsoid, compounded with the 
irrotational motion consequent on giving the shell additional angular velocities 

The velocity-potential of the irrotational motion will be 
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Hence, if h l9 h% 9 h 3 be the components of angular momentum, and M denote the 
mass of the fluid, p x its density, 
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If the system be disturbed from amotion of pure rotation, with angular velocity co, 
about the axis of % ; £, 77, I2 l5 2? fl 3 will all be small quantities, while £ will be 
approximately equal to co, and hence, on omitting small quantities of the second order 
and putting £ = <o in small terms, the equations of angular momentum, viz. : — 



h 



2 



ho — 



^ 3 p + /^r = 
h x q + ^ 2 r = 



r 



P 



}• where ^ ^ = O a + 77 

n 3 + c 



L 









MR. S. S. HOUGH ON THE OSCILLATIONS OF A 



become 



>2 __ *& . 



A (a, + .ft + v |^fl 1 + |M(/3* + y 2 K 



/S 3 + 7 s 



or — 7" 



B ( n a + *?) + /* a s + 7 2 



Xl 3 +|M(« s + r 3 )77 o) + (n 3 + r ? )[C + lM( a 3 + i8 3 )]^ = 



* * 

B (ft 3 ~4~ ^7/ H - * /^ 



a? — 7 3 



a 3 + 7 : 



2 



3 + 1M (y 3 + a 3 ) V 



/3 2 - y 3 



(fl 1 + 0[C + lM( a 3 +/8 3 )] w + A(n 1 + f) + v%--* i fl 1 +£M08HY)£ 



o" -4- <y 



CO 



= 



C(fl 3 + C) + IM^j^-'o, + iM(a* + ^) £"= 0. 
Helmholtz's equations of vortex motion are 



£-- 



2ofi 
a 8 + 7 a ft) °^ ^ 



2/3 3 



^2 _j_ fy2 



ct>ft l5 £ = 0, 



Hence, if we put g l9 tj 19 £l l9 ft 3? ft 3 , £ — a> each proportional to e ix \ and introduce, 
for brevity, the notation A! = -g-M (/3 3 + y 3 ), &c. 
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On expanding and arranging according to powers co, X, this determinant 
reduces to 

X* [A 08» + y 3 ) + v o* - y 3 )] [B («» + /) + /. («» - y 2 )] 

- o> 3 \ 2 [(a 3 + y 2 ) (£ 3 + y 2 ) (B - C - v) (A - C - p) + 4a 3 /3 3 (A + v) (B + /*) 

+ 4/*0Y (B - c - v) + 4wy (A - c - /»)] 

+ a, 4 {4« 3 /3 3 (B - C - i») (A - C - /*)} = (5). 

§ 2. Case of Shell without Inertia. 

If the shell be so thin that we may neglect its inertia compared with that of the 
fluid, we may put A = B = C = 0, and equation (5) then becomes 

(«» _ y 2) (£2 _ /) \ 4 - W 3 X 3 {/ - 3 (a 2 + /3 3 ) y 3 + 5a 3 £ 3 } + 4ft> 4 a 3 £ 3 = (6) ; 

when the system is symmetrical about the axis of rotation a? = /3 3 , and this equation 

reduces to 

X* (a 3 - y 3 ) 3 - ft> 3 X 3 (a 3 - y 3 ) (5a 3 - y 3 ) + 4ftAx 3 /3 3 = 0, 

the roots of which are 

„ a. 2 r, , //9««-7»\i8 

x3 = 4 l 1 ±V(^r 7 y 3 j} W- 

These are the same as the values obtained by Bryan ('Phil. Trans./ 1889, A, 
p. 208), for the case of a spheroid whose surface is free. As is there indicated, the 
modes of oscillation corresponding to these periods are such that the surface of the 
spheroid maintains its shape, but changes its position. Such oscillations will, 
of course, not be affected by supposing the fluid contained in a rigid shell without 
inertia, and we might have expected to obtain the same values for the periods, when 
the figure of the shell agrees with a possible figure of equilibrium of the fluid rotating 
freely. 

From (7) we see that the roots, if real, are positive ; in order that they may be 
real, we require that 9a 3 — y 3 and a 3 — y 3 must have the same sign. 

Hence a necessary condition for ordinary stability is 

y 3 > 9 a 3 or y 3 < a 3 , 

i.e., y must not lie between a and 3a. 

Returning to the case where a 3 ^fc /3 3 , in order that the roots may be real and 
positive, we must have 

(1) (y 3 - a s ) (y 3 ~ £ 3 ) > 0. 

(2) y 4 - 3 {a? + /3 3 ) y* 2 + 5a 3 /3 3 > 0. 

(3) {y 4 - 3 (a 3 + /3 3 ) y 3 + 5a 3 £ 2 } 3 - 16a 3 £ 3 (y 3 - a?) (y 3 - /3 3 ) > 0. 

MDCCCXCV. — A. 3 P 
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The 1st condition requires that y 3 should not lie between a 3 and /3 3 . 
Now 

y 4 _ 3 ( a 2 _j_ £2) y 3 + 5a *pi _- ( y 2 _ g a 2) ( y 3 _ 3/ g2) _ 4^8 

= (3a 3 - y 2 ) (3£ 3 - y 3 ) - 4a 3 /3 3 . 

The 1 st form shows that condition (2) is certainly satisfied if 

y 2 > 5 a 3 and also > 5/3 3 . 

The 2nd shows that it is satisfied if y 3 < a 3 and < /3 2 . 
Lastly, 
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Hence condition (3) will certainly be satisfied if y 2 > 9a 3 and also > 9/3 3 , or if 
y 2 < a? and also < /3 3 . 

Thus the roots of (6) will both be real and positive if 

y < a and also < /?, 
or it 

y > 3a and also > 3/3. 

These conditions are sufficient, but not necessary, to ensure stability ; the neces- 
sary conditions are given by the inequalities (1), (2), (3). 

The analytical conditions here discussed are approximately realised in the case of a 
liquid gyrostat (vide ' Nature/ vol. 15, p. 297) mounted on gimbals in such a way that 
the centre of gravity is held at rest. The inertia of the gimbal-rings will be unim- 
portant when the rotation is rapid, and, if we may also neglect the inertia of the case 
compared with that of the fluid, the gyrostat will be stable when set rotating about 
its least axis ; it will also be stable when set rotating about its greatest axis when 
this axis is, at least, three times as great as either of the others. It will, however, 
certainly be unstable when set rotating about its mean axis. 

§ 3. Approximate Solution of the Period Equation. 

Let us for the future suppose that the cavity which contains fluid is approximately 
spherical, so that — — ? — are small quantities, 

1 ry ry L 
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Suppose 

a — cy ft — ry 

~ — €1, z:z 69. 

y i y 

From (3) we have 

therefore, if we neglect squares of e l9 e 3 , 

fV € i = v\H = ?C say. 

In the case where the thickness of the shell is finite compared with its linear 
dimensions q will be a finite quantity; when the shell is thin q will be large, and 
when the fluid nucleus is small compared with the dimensions of the whole system q 
will be small, the densities of the fluid and of the crust being supposed comparable 
with one another. In all cases q will be positive. 

If e l9 € 2 are each equal to zero, the equation (5) becomes 

ABX 4 - [(C - A) (C - B) + AB] W + (C - A)(C ~ B) o 4 = 0, 
or 

(X 3 - *,*) [ABX 3 - (C - A) (C - B) a> 3 ] = 0. 

Thus we obtain as a first approximation to the roots 
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Next let us retain first powers of e 1? c 2 in (5) ; this equation then becomes 
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+ € 2 [ABX 4 - *> 3 X 3 {(0 - A) (C - B) + 2AB + qBC] 

+ co 4 (C~ A){2(C-B) + 5fC}] = . (8); 

dividing by ABX 3 — (C — A) (C — B) o> 3 , and putting X 3 = co 3 in the terms which 
contain e x or e 3 as a factor, we obtain as a closer approximation to the root X 3 = co 2 , 
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Again, dividing by X 3 — co 3 in (8), and putting X 2 = — — — - ^ a> 3 in the small 

terms, an approximate value of the second root, correct to first powers of e l9 e 2 , is 

given by 
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to the same order of approximation. 

This approximation involves the assumption that e v e 2 are small compared with 

T-— — , —Z — ; the approximate value of the root will, however, be the same if we 

q A C — B 

suppose — - — , — ^— to be small quantities of the same order as e l9 e 2 . 
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Retaining only finite terms in (5), we obtain as a first approximation to the roots 
X 2 = o) 2 and X 3 = ; also, the independent term in (5) is a small quantity of the 
second order in k 19 k 2 , e 19 e 2 . Thus the root which approximates to X 3 = will be of 
the second order. Regarding X 3 as of the second order, and retaining only terms of 
this order in (5), we get 



X 3 . [4ABa 3 /3 3 ] - a) 3 . 4a 3 £ 3 . {Bk 2 + qGe^} {Ak x + qCeJ = 0, 
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or 

C 



X s = a>*. (^ + gr— €l j (% + 2 B eg 

= o> 2 . (k x + g^) (k 2 + ge 3 ) 
to the same order, and this is the value obtained above (10). 

§ 4. Application to the Case of the Earth. 

The nature of the two types of oscillation will be found fully discussed in the 
Appendix. It is there shown that the oscillation corresponding to the root o> (1 + E) 
is that previously examined by Hopkins in his ' Researches in Physical Geology/ 
whereas the second type is analogous to the motion of a rigid body when disturbed 
from a motion of rotation about a principal axis. 

If the Earth could be regarded as a system such as we have been considering, we 
see that in addition to the ordinary Solar and Lunar Nutations, which would be of 
the same nature as when the Earth is supposed solid throughout, there might exist 
certain free nutations the amplitude of which could only be determined by observa- 
tion. If the amplitudes were sufficiently large, the oscillations corresponding to the 
root X = o) (1 + E) would render themselves visible in the same way as the Solar and 
Lunar Nutations, namely, by small periodic displacements common to all stars. The 
period of these displacements would be 1/E sidereal days, and a knowledge of it would 
enable us to determine E, a quantity which depends on the form of the internal 
surface and the thiokness of the crust. 

The oscillations which correspond to the root X = ^-v/U^i + qe^ (k 3 + qe 2 )} 
would manifest themselves in a different manner. They are, in fact, similar to the 
"Eulerian" nutation (vide Tisserand, ' Mecanique Celeste/ vol. 2, p. 494), and will 
involve a small periodic change in the latitude of places on the Earth's surface, as 
found by meridian observations of a circumpolar star, this change taking place in a 
period of {(k x + qe^ (k % + ^ 2 )}"~ l sidereal days. 

Now it appears probable that in oscillations of long period, such as Precession, the 
effects of fluid friction would be to make the internal fluid move with the crust 
as if rigidly connected to it (Tisserand, ' Mecanique Celeste/ vol. 2, p. 480, or 
Lord Kelvin, 'Popular Lectures and Addresses/ vol. 2, p. 244). Hence, if &, 
^J, & be the principal moments of inertia for the Earth as a whole, supposed 
symmetrical about its axis of rotation, the Theory of Precession will still enable us 



to determine the value of — ^ — , as 3-^5. 



But, if we put ac 3 = k 19 € 3 = e l9 and denote by M the mass of the fluid, 



<& = C + f Ma* = C {1 + q (1 + 2e x )}, 

i = A + |-M(a 3 + r 3 ) = A + qC (1 + € 2 ). 
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<jp ~ g __ C — A 4- qGe^ _ k x -f ^ 
CD — 0(1 +7+2^) ~~~ 1 -f gr 



very nearly. 

Therefore the period in which the latitude variations will take place is 



sidereal days = 



*i + <Z 6 i 



©-&. 1 + q 



sidereal days ; 



when the Earth is supposed solid throughout, this period is 



© 



GC~® 



sidereal days. 



We thus see that if the Earth consisted of a rigid shell containing a homogeneous 
fluid nucleus, the theoretical period of 305 days, calculated on the assumption of the 
Earth's rigidity throughout, would be diminished in the ratio 1 : 1 -f g, where q is 
an essentially positive quantity, whose magnitude increases with the size of the 
nucleus. 

In order to form some idea of the magnitude of this effect, let us suppose that the 
fluid and the crust have the same density p, and that r, r 3 are the mean radii of the 
fluid nucleus, and of the Earth as a whole. 

We then have approximately 



li = y = -£ 5 - 7rpr b € l and C = yg irp (r x B — r 5 ). 
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and the period will be diminished by ( — j X 305 days. 

Taking the mean radius of the Earth as 4000 miles, we obtain the following table, 
where the .first line gives the thickness of the crust in miles, and the second the 
diminution of the period in sidereal days : — 



Thickness of crust in miles . . -. . 


2000 
5 


1000 
21 


500 


250 
221 


100 


Diminution of period in days . . 


156 


269 
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Now, although in the light of Professor Newcomb's work (' Astron. Soc. Monthly 
Notices/ March, 1892), it appears probable that this effect would be modified by the 
elasticity of the Crust, it could scarcely be reversed if the fluid nucleus were of any 
considerable extent. We must, therefore, conclude that the observations on latitude- 
variation, so far from establishing the existence of a fluid interior, as supposed by 
M. Folie, rather tend to confirm the views hitherto maintained by physicists on 
other grounds, that there can be no internal fluid mass of any considerable extent. 



APPENDIX. 

Treatment of the Problem by Lame Analysis. 
§ 1. Equations of Motion of Fluid. 

Let us refer to rectangular axes rotating with angular velocity o) about the 
axis of z. The fluid is supposed to have no motion relatively to these axes other 
than that due to the small oscillations with which we are dealing. 

Let u 9 v, w, be the velocity-components at any point x, y, z relatively to these 
axes : we shall, as is usual in small-oscillation problems, neglect squares and products 
of the small quantities u, v, w. 

The actual velocity-components parallel to the instantaneous positions of the moving 

axes will be 

u — a>y, v + o)X > iv, 

and the differential equations of motion of the fluid are therefore (Basset, i Hydro- 
dynamics/ p. 22) 

- -a>(v+a>x)-o>v = ~- dx (Y, - _ 



dt dz \ p } 

where V x is the gravitation-potential of the forces to which the fluid is subject, p the 

fluid-pressure, and p } the density. 

Putting 

xjj = Y l -p/ Pl + !>(*>* (x* + y*) ....... (1), 

the above equations reduce to 

du/dt — 2o)V = dxfj/dx"" 

dv/dt + 2am = Sift/dy }>.... (2). 

dw/dt = dxfj/ dz 



480 



MR. S. S. HOUGH ON THE OSCILLATIONS OF A 



We have, in addition, the equation of continuity 



du/dx + dvfdy -f div/dz = (3). 



Equations (2), (3) are sufficient to determine u, v, iv, xfj, subject to certain boundary 
conditions. 



From (2) we obtain 
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Applying the operators 3/3<£, 3/3;?/, 3/3^ and adding, we obtain by means of (3) 
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3 3 3 2 3 3 ' 
- 3^ 3 T 3?/ 3 T 3* 3 



or, by the third of equations (2), 



4W & 8 ^ 3f 2 l w 



= . 



(5). 



This is Poincar^'s differential equation for the oscillations of a mass of fluid about 
a steady motion of pure rotation. 

Let us now suppose that the system is executing one of its component harmonic 
vibrations. 

Assume that 



and 



u == u Y <£ x \ 



v = v ± e ix \ 



\jj = xjj^K 



w 



w x e 



i\t 



Putting these values in (4), (5), and dividing out by the time factor, we get 



u 



v, 



V)y = 



4o> 3 ™ X 3 
J. 

4*> 3 .- X 3 
1 3ifo 



•x 3 ^3 



3 fl 



£A TT - "T- 2&) ~ 

OA' ' oy 






2oj 



dx 



r 



» » « 



• • ( 6 )> 



^X 3; 



.2 



while i|r T satisfies the equation 



3 2 ^i , 3 3 ^ 
3^ 3 3?/ 3 



+ i 



4 w 3 \ 3^i 
V 7 df 



= 



e • 



• • 



• (7). 
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§ 2. The Boundary Conditions, 

The position of the shell at any instant may be defined by means of three coordi- 
nates, l9 2 , 3 , which denote the small angular displacements, about the axes of 
reference, of the shell from the position it would occupy in the steady motion. 

The displacements parallel to the coordinate axes of the point of the shell, whose 
coordinates are (x, y } z), are 

— y0 s + z0 2 , — z0 x + &0 8 , — x0. z + y0 v 

If cos a, cos/3, cosy be the direction cosines of the normal to the undisturbed 
surface, the normal distance between this surface and the displaced surface will be 

(— y0 s + 20 2 ) cos a + (— z0 x + x0 s ) cos ft + ( — x0 2 + yOi) cos y 

= 9 X (y cos y ' — z cos ft) + #2 (% cos a — x cos y) + # 3 (x cos /3 — y cos a). 

The condition to be satisfied at the boundary is that the rate of increase of this 
length must be equal to the component velocity of the fluid, relative to the moving 
axes, in the direction of the normal to the undisturbed surface. Now as these rela- 
tive velocities are all small quantities whose squares we are neglecting, it is 
unnecessary to distinguish between the velocities at the disturbed and undisturbed 
surfaces ; thus, at the latter surface we require 

u cos a + v cos /3 + w cos y — x (y cos y — z cos /3) 

+ 0% (z cos a — x cos y) + #3 (# cos fi — y cos a). 

Putting X — 0\e ixt , &c, and omitting the exponential factor, we obtain 

u x cos a + v Y cos ft + Wi cos y = iX [0^ (?/ cos y — s cos /3) 

+ &% (z cos a — x cos y) + #'3 (# cos /? — y cos a )] ; 

or, putting in the values of w x , ^ 1? w x , ^ from (6), 

7~i ^ I -^- cos a + -^ cos p + -3" cos y 1 — I 

4cm z — V [ ax vy oz l \ A, 3 / 

2coi 



4a) 2 - X 3 



cos a ;— cos B 



. dy 

= X [#'i (2/ cos y — 2 cos /3) + 0' g (z cos a — a? cos y) 

+ 0' 3 (a? cos fi — <y cos a)~] . (8). 

Let us now put 



X 3 

MDCCCXCV. — A. 3 Q 



4w 2 3 , 
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Corresponding to any series of points whose coordinates are denoted by (x, y, z), 
we shall obtain a new series whose coordinates are (x, y, z), which will be real or 
imaginary according as X s is greater or less than 4oA We will take as the standard 
case that in which \ 2 > 4o) 2 . 

If the point (x, y, z) trace out any surface whose equation is f(x, y, z) = 0, the 
corresponding pointy, ;?/, 2') will trace out a new surface whose equation i$f(x } y,rz) = 0. 
The part of this latter surface which corresponds to the real part of the surface 
f (x, y, z) = will, however, be purely imaginary if X 2 < 4a) 3 . 

If (cos a, cos /?, cos y), (cos a, cos /3', cos y) be the direction-cosines of the normals 
to the two surfaces, we have 

cos a : cos /3' : cos y = dfj dx : dfj dy : dfj dz' 

— dfj dx : dfj dy : r (dfj dz) 
= cos ol : cos /3 ; r cos y. 

Substituting in equations (7), (8), the differential equation for fa takes the form 

d*fa/da? + Wfa/dy* + d*fa/W* = ...... . (9), 

while if / (x, y ,z) = be the equation to the undisturbed boundary of the fluid, ^ 
must satisfy the equation 



. i 



dfa f dfa , dfa 

cos a + -f cos B + ^ cos y 



- — 2coi < 



(ico* - X 2 ) \ . 



*l' 



# 



r 2 ' I T2; ' cos a - 



dfa t 

T cos a 



. ty 



cos y' — tz' cos /3' 



/ 



<AJ f 

— cos y 



-f- 3 ' (x cos /3' — y cos a') __ 



dfa 



IX 



cos /3' 



* # • 



. . (10) 



at the surface f(x, y, z) = 0. 

The problem of finding the motion of the fluid is thus reduced to that of obtaining 
solutions of equation (9) consistent with the boundary condition (10) at the surface 
f(x, y, rz f ) = 0. 

§ 3. Case of Ellipsoidal Surface. 

Hitherto, no assumption has been made as to the form of the surface of the fluid. 
Let us now suppose that it is given by the equation 



X" 



, if 



p % + p 3 _ I 



O 4* 



Hi 



e 3 



-L • 



. . . (A), 
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that is to say, that it is an ellipsoid whose principal axes coincide with the axes of 
reference. 

We will take as the standard case, that in which 

p 2 > c 2 > ¥ > 0. 
Put 

p % — <? = (p S — C' 3 ) T 3 . 

The equation to the auxiliary surface f(x, y, rz f ) = becomes 



.,2 „& ~'2 



_. I _1____ J I -1 /A'\ 

p 2 ' p 2 __ £2 ' p 2 __ c '2 X • • • \ A /• 

Let us now introduce two sets of elliptic coordinates (p 9 p, v), (p. //, v), connected 
with (x } y, z) (x, y, z') respectively by the equations 

x fiv y \/{i $ — • & 2 ) (W - — • p 3 )" z \/ {f — • yu, 2 ) (c 2 — v 2 ) , * 

~~p~~U> vV^j ~ 67(7^ 6») ' ^/(p 2 - c 2 ) ~~ <? v /(c 2 - Z> 2 ) * ^ ^ 









p' will be equal to p for points which lie on the surfaces (A) (A'), but not otherwise. 
Let us also put 

X = x/p, Y = <//vV ~ n Z = zlJtf - c») = z'ly/{p'* - c'% . (13) 
for points on these surfaces ; so that X, Y, Z are subject to the relation 

X 2 + Y 3 + Z 3 = 1 . (B). 

X, Y, Z may therefore be regarded as the coordinates of a point lying on a sphere 
of unit radius. 

Denote by R, M, N three conjugate Lame functions of the elliptic coordinates 
p, fi, v, and by R', W, N' three similar functions of the coordinates p\ p! } v. 

A form of solution of equation (9) convenient for satisfying boundary conditions at 
the surface (A') is 

i/r 1 = sA'R'M'N' (14). 

The effect of the fluid on the motion of the shell will depend only on the fluid 
pressure over the surface, and this by (1) will involve the value of \\s x at the surface. 

3 Q 2 
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To find the value of \}j } at the surface (A), we may transform the expression (14) for 
i// x first to the surface (B) and then from (B) to (A). 

Now, by a known property of Lame products (vide Heine, ' Kugelfunctionen/ vol. 1, 
§ 89), if M, N be two conjugate functions of order n, the product MN at the surface 
(A) will transform into a surface harmonic of order n at the surface (B) ; and, con- 
versely, any surface harmonic of order n at the surface (B), when transformed to the 
surface (A), can be expanded in a series consisting of Lame products with constant 
coefficients, each of which products will be of the n th order. 

The same conclusions will hold for the surface (A') and the sphere (B). 

We can thus express the value of \p x at the surface (A) in terms of a series of 
Lame products, in which each term will be of the same order as that from which it 
arises in (14). 

The couples on the shell due to fluid pressure are 



p dcr (y cos y — % cos j3), \p dcr (z cos y — x cos /3), pdcr(x cos /3 — y cos y) 



where dcr is an element of the surface and the integrals are taken over the whole 
surface. 

If P denote the perpendicular from the centre on the tangent plane to the ellipsoid 

P 

(A) and I = —^—==~ 



y cos y — z cos ft = Yyz 



IP' 



&■ 



3 



¥ 






Vyz (c 3 - ¥) 

(p2 ^( p 2^ 



and Vyz is proportional to 



y^ _ tf) \c* _ p*) y[p _ ^) ( c * ~~^) - ZM X N 1? 



where M 1? Nj are two conjugate Lame functions of the second order. But, if MN be 
any two conjugate Lame functions different from M 1 N 1 , U^MNM^N^ dcr = 0. For, if 
we transform to the surface of the sphere (B), I dor is equal to the corresponding 
element of the spherical surface, and MN, M X N X transform into two different surface 
harmonics. 

Thus the only term in \jj 1 which can give rise to any couple about the axis of x will 
be the term involving the Lame product M^N"-^ 

Similarly the terms which can give rise to couples about the other axes will be of 
the second order. These, as we have seen above, all arise from terms of the second 
order in (14), and, in order to evaluate these couples, it will be unnecessary for us to 
calculate any coefficients in \\s l other than those of terms of the second order. 



ROTATING ELLIPSOIDAL SHELL CONTAINING FLUID. 485 

§ 4. Transformation of Boundary Equations. 

Let us now transform our boundary conditions to the surface of the sphere (B). 
We can at once express the right-hand member of equation (10) in terms of 

X, Y, Z ; for 

, P'a? P'X 

cos a = 

COS ft = 

cos y = 

where P' has the same signification with reference to the ellipsoid (A') as P has with 
reference to the ellipsoid (A). 

The right-hand member therefore becomes 



P 3 

Yy 


~ P ' 




" V(P* - & 2 ) ' 
P'Z 


p 2 - c' 2 


~ vV - ^) ' 



- XV . P' . 



r% — 7)3 f — /»2\ 7,2 
nr ° ^ VyT, /?' V _ / 7V i /Q> YV 



(15). 



Consider next a single term of order n in \fs l9 say i/^ = R'M'N'. We have seen 
that M'N' is expressible in the form cS„, where e is some quantity which does not 
vary over the surface of the sphere (B) and S* is a spherical harmonic function of 
degree n in X, Y, Z. 

If dn r denote an element of the normal to the surface ( A'), we have P' dn = p dp\ 
and therefore, 

-^ COS a + — COS p -*- -^7 COS y = "^-7 = — , 777 . 
d# oy as / o?i p d/> 



Now M', N' are independent of p. Therefore, w T hen \jj x = R'M'N 



if- cos a + ■— - 1 cos p 4- -77 cos y = -7 rrlr) = — a 
dx oy os ' p op p op' 



Next let ds be an element of a line through (x, y, z) whose direction-cosines are 

COS Q COS OL \ 

_ __ — £L - r — ) , and which, therefore, lies in the surface (A'), 
sin 7 sm 7 / 

Then 

-— cos <x ~ cos p =smy ^f- , 

d?/ da? ' ds 

and when ifo = R'M'N' = R'eS tt , 

"3s 1 == '& ^ ^ * ^ "^ ^ e * "sf * * * * * * * ^ ^* 
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Now, the element ds lies entirely in the surface (A'), and corresponding to every 
point of it there will be a point X, Y, Z which lies on the sphere (B). Hence, as 
B/, e do not vary over the surface of the sphere, they will remain constant as we 
pass along the element ds. Thus we have 



£ < R ' e > = ° ; 



and, therefore, from (17) 



Cs 



'^tl — TV ^ 

as as 



K e^ _ . ls + dY . 3g + 3Z . ^ 



[3X' 3^ 



n 



But from (13) we see that 



1 a 



w 



p ds 



dy 



cX 

ds 

3Y _ 

ds 

az 



cos /3 /N 
sin 7' 



PV 



p (p 2 - & 3 ) 



1 /cos a r 



cosec y 



P'aj 



P'Y 



/V(p 2 -& 3 ) 



cosec y 



pV(p 8 -i b ) 



cosec y = 



FX 

/> \/0> 3 - & 3 ) 



cosec y 



3; 



Z 



= 0. 



Therefore, 



and 



3^ 
3y 



cos a 



9^i 

3,9 



df } 



IV 



F.K'e 



dS^ 1 



— x <xtt — J x — — Y — i cosec V 



cosjS 



'_ • '?&_ F - R/ 



Smy 35 _ p^/(p»-J*) 



— < 



dY ~ dX 



■ ( 18 )- 



From (15), (16), (18) we see that the boundary equation, on the assumption of the 
form (14) for n\t v takes the form 



\. 



dp 

- - \ z Ae\ 



2m\tk! ~^ ¥ (X^ 



Y 



n 



ax 



> 



p (c a - P) 
vV - V) (p* - c 2 ) 



Y^j "T" P 2 — tt 5 rr ZX -J- 6 



W 



s/(p* ~ c>) 



3 . /(r» 



vV - ¥) 



XY 



(19). 



Now the function X ^~ — - Y ^=p is itself a spherical harmonic function of order n 9 



as, 



and both S* and X^~ — Y^£ may be expressed as linear functions of the 2n -f 1 
independent harmonics of the nth order. 
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Omitting for the present the terms of the second order, if we equate to zero the 
coefficients of the 2^ + 1 independent harmonics, we shall obtain a series of 2n + 1 
linear equations connecting the 2n + 1 quantities A' which occur with Lame products 
of order n. These equations show that all the quantities A' vanish except for certain 
values of X which satisfy the determinantal equation obtained by eliminating them. 
The roots of this equation will determine the possible periods of free oscillation, and 
when the system is oscillating in the mode corresponding to any one of these roots, 
Lame products of one order (n) only will appear in i^. 

These modes of oscillation do not involve any motion of the shell, and it is evident 
that they could not be generated or destroyed by any disturbance communicated to 
the shell, if the fluid be free from viscosity. 

We proceed now to examine more closely the modes which depend on terms of the 
second order. As we have seen above, terms of different orders correspond to different 
fundamental modes ; and therefore we may for the future suppose that the second 
order terms alone exist in xjj^ 



§ 5, Lame Functions of Second Order. 
A Lame function of order n is a function R of one of the four forms 



where V n denotes a rational, integral, algebraic function of p of degree n, and R 
satisfies the differential equation 

0* - b>) (? - c*)^ + p(2 P * -V-c*)f p = [»(n + 1) p 3 - B] E . (20), 

where B is a suitably chosen constant. 

The Lame functions of the second order are, therefore, the three functions 



px /~pF^¥, ps/fF-c\ </\p* - W) {? - <?) . . . (21), 

together with two functions of the form p % + /3 (22). 

To find these latter functions, substitute in equation (20) with n = 2 ; we obtain 

2 (p 3 - Z> 3 ) (p 3 - c 3 ) + 2p 3 . (2p 3 -¥- c 3 ) = (6p 3 - B) (p 3 + /3). 
Equating coefficients of p 3 , and the terms independent of p, we get 

- 4 (6 3 + c 3 ) = 6/3 - B, 26 3 c 3 = - /3B. 
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Eliminating B, the values of ft are given by 



3/3* + 2 (6 3 + c 3 ) /3 + 6V = . 



5 6 $ * * 



. (23). 



Let us now apply the formulae (16), (18), to the different forms of B' ; take first 

At the surface 

M'JN' = j»V (/*' a - & ' 2 ) • *V (^ 3 - ^ 3 ) = ^JTTTm x y = ^V (c' a - & a ) • XY. 



Therefore, when ft = B'M'N', 
from (16) 

-^- cos a + -H 1 cos « + ^r cos y — r . — fr-r — — 6^c */(c 

3a; 3y ^ 3z' ' p^ (p 3 — 5 3 ) v v 

from (18) 



■2 



¥) . XY 



"^ 



y (24). 



^ cos a! - ^- l cos /3' = F 6V*/(c' 3 -<&*) . (X ! 

dy da? v x 7 x 



2 



V2\ 



J 



Similarly, when B' = P V (p* - c*), 



M'N' = 6c' 3 y (c' 3 — 6 3 ) XZ at the surface. 



an 



d 



-~ — COS a -J — ^— cos p + -w-f cos y 

OtX/ 02/ oz 



T)/ 



2/) 3 







'2 






COS a - — ^ cos p 






jv &c'* v 7 (^ - & 3 ) XZ, 



7^? 



F Vjnrp(- YZ ) 



\-«o) 5 



and when R' = \/(p' 3 — - & 3 ) (// 3 — c' 3 ), 



M'N' = 6c' (c' 3 - 6 3 ) YZ, 



-7T-* cos a 4- -Tp- cos B + -rr COS y 

ox> d'u dz ' 



P' 



9/»2 — 7»2-_i. c '3 



2p 3 -5 



3y 



V^"( p 8 - 5 3 ) (p 3 - c*) 



:/r & c ' (c' 3 - & 3 ) YZ, 



cos a 



3$ 



A /Y/,2 _ JZ\ 

cos /3' = P' ^ °J he' (c' 3 - 6 3 ) XZ 



>*■ (<4u/« 



Take, now, the form (22). 
When 



E = (p 3 + /3), RMN = (p 3 + /3) (/x 3 + /3) (*> 3 + /3) ; 
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but p 2 , p?, v* are the three roots of the equation in a 






therefore, 



(a — p 2 ) (a — fi?) (a — -*i> 2 ) = a (a — 5 2 ) (a — c 2 ) — cc 2 (a — 6 2 ) (a — c 2 ) 

— y 2 ct (a — c 2 ) — 2 2 a (a — 6 2 ) ; 
putting a = — f3 in this identity, we obtain 

RMN = (/3 + p 2 ) (/3 + />t 2 ) (j8 + *> 2 ) = £ (0 + W) (f3 + c 2 ) 

+ x 2 (f3 + 6 2 ) (/3 + c 2 ) + ^ (£ + c 2 ) + 2 2 /3 (/3 + 6 2 ), 

and, therefore, at the surface 

RMN = (X 2 + Y 2 + Z 2 ) /3(/3 + & 2 ) (/3 + c 2 ) + p 2 X 2 .(/3 + & 2 ) (£ + c 2 ) 

+ (/> 2 - 6 2 ) Y 2 . /3 (£ + c 2 ) + (p 2 - c 2 ) Z 2 ./3 (/3 + 6 2 ) 
= (/>*+ /3) [(/3 + & 2 ) (/3 + ^^^ 

or, since by (23) (/3 + 6 2 ) (J3 + c 2 ) = — /3 (y8 + c 3 ) - (3 (p + 6 2 ), 

RMN = (p 2 + )8) [fi (0 + c 2 ) (Y 2 - X 2 ) + (3 (f3 + ¥) (Z 2 - X 2 )], 

therefore, when R' = p f% + /3', 

M'N' = p (P + c' 2 ) (Y 2 - X 2 ) + p ((3 f + ¥) (Z 2 - X 2 ), . 



and if ^ = R'M'N' 

-~r cos a + -^f" cos p + -^r cos y 

doc ay cte ' 

= 2F [# (/3' + c' 2 ) (Y 2 - X 2 ) + £' (& + ¥) (Z 2 - X 2 )] . (27), 
and 

d ~P cos a' - ^ cos 0' = F £ + /m 2XY {2/3' (£' + c' 2 ) + P (P + & 2 )}, 
or since (3' satisfies an equation similar to (23), 

^ft cos *' '- gfr cos ff = - F ■ ^ + ^ - 2XY {(c' 2 + P) 5 2 } . . (28). 
Let us denote by /3' l3 /3' 2 , the two values of P, and assume that 

MDCCCXCV. — A, 3 R 
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^ = A, (p' 2 + fi\)(^ + PW* + ? x ) + A 2 (p' 2 + /3' 2 )(p' 2 + P\)(v'> + #,) 






+ 






Bi 



-Do 



; p' vV 3 - 6 s )./*' vV* - &V' y(6 a - ^) 



p' vV 2 ~ c'V x/(c' 3 - p'V ,/(c'*- i/*) 



J. c ' ( c '3 _ J9) 



- ^/(p'a _ 6 2 )(p' 3 - c' 2 ) V(p' 2 - 6 2 ) (c' 2 - p' 2 ). \/(o 2 - */ 2 )c' 2 - i/ 2 ). 



A. 



From equations (24) . . . (28), we see that (19) now takes the form 

2 A, {j3\ (^ + c' 2 ) (Y 2 - X 2 ) + /3\ (ft 4- 6 2 ) (Z 2 - X 2 )} 
+ 2A 2 {ft (ft + c' 2 ) (Y 2 - X 2 ) + ft (ft + & 2 ) (Z 2 - X 3 )} 

2p2 _ J3 _ </2 



i T> _jg_ ^L yy «l "R -"£ !—■- YZ 4- "R 



\//3 3 — & 



p v//3 a — C 



f"2 



\/(p 2 - 6 3 ) (p 2 - c' 3 ) 



— YZ 



2 cm 



P s/(P % - V) 



- 2A 1 6 2 (^' 1 + c'%^ + (B\)XY - 2 A g 5 2 (ff a + c' 2 )(p 2 + /3' 2 ) XY 
+ BjpN/p^^P (X 2 - Y 2 ) - B 2 pN/p 2 - c' 2 YZ 

+ B 3 y/Q^8*)fiF=c*y. XZ 



- XV. 



^i 



^/(pS _ &2) (p« _ cS) 



YZ + 0' 



— c 



2 



P \/j0 3 



zx + 0'< 



5 3 



P vV — ^ 3 



A. X 



(29). 



6. Calculation of Coefficients in \jj lt 



Equating coefficients of Y 2 — X 3 , Z 2 — X* 2 in the two members of (29) 



(Dl 



KP\ (Pi + c' 2 ) + A 3 /3' 3 (#, + c' 2 ) + ^ Bi = 



> . 



» • 



A^ (/3\ + Z> 2 ) + A# 2 08', + 6 2 ) 







. (30). 



Multiply the first of these by 2 and add to the second. Reducing by means of 
(23) we obtain 



°ai 



b\ [- A, (c' 2 + p x ) - A 2 (c' 2 + fa)] + ^ B, = .... (31). 



P 

Multiply by -—- and subtract from the first of equations (30) ; then 



0)1 



A (p 2 + /3'0 (c' 2 + /3',) + A 2 (p 2 + fa) (c' 2 + /3' 2 ) + y f 1 



2p^ 

~^/ j:> i 



B, = 0. 
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But if we equate coefficients of XY in (29) we obtain 



4g)z 



B x (2p 2 - V) + ~b\ [A x ( P 2 + PJ («* + fi\) + A 2 (p» + ? % ) (c' 2 + /3' 2 )] = - X 2 rW 3 . 



Hence 



or 



or 



4o) 2 



Bj (2p 2 - fc 2 ) + — Bj (6 a - 2/a 3 ) = - \W 3 Z> 2 , 



X 3 



Bj (2p 2 - ¥) t 2 = - \W 3 6 2 , 



B 1= = - 



\ 2 & 3 



V - & 



0'. 



3 " 3 



(32). 



Again, equating coefficients of XZ, YZ in (29), we have 






2 



c 



/>\/V - c3 ) 



2p » - y - g on yv - ^) _^ 

3 \/(/> 3 - fc 3 ) (p 3 - c' 3 ) T * k ^ ^a __ y) — 



/•2 __ A3 

_ XV #' - — — — ■ 

1 vv - & 3 ) (p 3 - c 3 ) ' 



or, since p 3 — c 3 = r 3 (p 3 — e' 2 ) 



B 3 (2/> 3 - c' 3 ) - 



2on 



B 3 (p 2 - c' 2 ) = \ 2 T0' a c 2 



9 



B 3 (2p 3 - & 3 - c' 3 ) + ^ B 2 (p 3 - c/ 3 ) = - X 3 t^! (c 3 - ¥) 



X 



.J 



. (33). 



Now, as we have seen above, 

Ai 0>' 2 + ^) 0*' 2 + £'i) („* + fiT,) + A, ( P ' 2 + #,) (/.' 2 + £' 2 ) (^ + #,) 
= A 1 {/3\(/r i +& 2 )(/r i + c ' 2 ) + a; 2 (^^ 
+ A 2 {/3' 2 (/3' 3 +& 2 ) 08VfC»)+s» (£' 3 +& 2 ) ( / 8' 3+ c' 2 )+3/ s /3' 2 (£' 3 +c' 2 )+s' 2 0' 2 (/3' 3 +& 2 )} 

which by (30) is equal to 



^iPMP\ + & 3 )+ A 2 ^' 2 2 (/3' 2 + 6 2 ),+ aV 2 { A 1 (fi\ + & 2 ) + A 3 (/3' 3 + ¥)} +y 



2 



Wi 



*i 



A. 

= A, (fi\ + V) { - f (& 2 + c' 2 ) £' x - Wo"} + A 2 (/3' 2 + & 2 ) { - | (& 2 + c' 2 ) £' 3 - |&V 2 } 



0)1 



+ — B x (a; 2 - y»), 



3 R 2 
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and this by means of (30) and (31) reduces to 



Hence the complete value of \jj l is 

— B, { - W + ^ ~ y*} + Bjxy + BX + B s yz' 



Ml 



Bi { — i& 2 + ^ 2 — y 2 } H- B^ + - xz + — ?/z 



. . (34) 



where B 1? B 2 , B 3 are given in terms of & v 0' 2 , 6\ by equations (32), (33). 

From equations (6), (34) it will be seen that, in the motions with which we are 
dealing, u, v, w are linear functions of x, y, z. Hence the components of molecular 
rotation of the fluid, which involve first differential coefficients of u, v, w, will be inde- 
pendent of x. y, z. This justifies the assumption made in § 1 of the paper. 

§ 7. Calculation of Couples on the Shell due to Fluid Pressure. 

At any point of the fluid the pressure is given by (1) ; we have, viz. :— 

P = Pi{Vi + W(*P + tf)} ~ Pit- 

Let us now refer to a new set of rectangular axes, Ox l9 Oy l9 Oz l9 coincident with 
the principal axes of the ellipsoid in its displaced position. The direction-cosines of 
one set of axes relatively to the other are given by the scheme 



Hence 





x l 


V\ 


■h 


X 


1 


-h 


<V 


y 


#3 


1 


-A 


z 


-*» 


o, 


1 



x — x Y — y x d z + z x B % , 

y — V\ — hPi + x \ e & 

z = z l —x^ + y x e v 



(35). 
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and, neglecting squares of V 2i 9 % , we obtain 

P = Pi [Vi + W fai* + Vi) - G>*y&A + «>%zAl - pi$ • • • (36), 

where in the small tQrm \\$ we may replace x, y, z by x l9 y v z lt 

If L, M, N denote the couples on the shell about the axes Ox l9 Oy i9 Oz l 



M = 



N = 



p da 
\p dcr 



tyih *Vi 



c A 



tf 



p 3 - c 3 



-6 



2 






,2 _ /.2 



(P 2 - W) (p 2 - C 2 ) 
-~c 2 [C 



] Ppy^i ^ 



p-(p 



pda< 



1 3^^ 1 tf^ 



H 7 



9 



2 /„ 2 



p 2 (p 



rr 

— "ft 3 ) J ^ > ^ 1 ^ 1 ^ 



where <icr is an element of the surface of the displaced ellipsoid, and the integrals are 
taken over the whole surface. 

Let us now consider separately the parts of these couples introduced by the 
different terms in the expression (36) for p. 

(a) Take p = /^V^ 

The pressure at every point is the same as if the fluid were at rest under a 
potential Y v 

V x will, in general, consist of three parts due respectively to (a) the attraction of the 
shell ; (/3) the mutual attraction of the fluid particles ; (y) any external attracting 
system. 

If the part (a) gave rise to any couple, it would be exactly counterbalanced by the 
couple on the shell due to the attraction of the fluid, since the attractions of the shell 
on the fluid and of the fluid on the shell are equal and opposite. 

The system of forces (/3) also form a system in equilibrium, and, therefore, can 
give rise to no resultant couple on the shell. Thus no couple can arise from the 
mutual attractions of the parts of the system. 

The pressure at any point due to the part (y) is the same as if the fluid were at 
rest. Thus the couples due to any external attracting system will be the same as if 
the fluid were supposed to be solidified. If we add to these couples, due to the 
attraction of the external system on the fluid, the parts due to the direct attraction 
on the shell, we see that the total couples due to any external system will be the 
same as if our system were solid throughout. 

(6.) Take 

p = Wpi i x i + Vi) - *>*pi*i (vA - a* A). 

Integrating over the surface of the ellipsoid 



x- 






V\ 



£2 



vrnmLmm 



z- 



2 



/3 3 — C 3 
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we have 



Therefore, 
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j JP do- y x V = A */> (p 2 &2 ) f (p 2 - ^ 2 )% &c. 
J jpP rfcr 2/i^i = - <» 2 0iPi • A «"/> (p 3 - &f (p* - c 3 )% 



and the corresponding part of L is 

- TS TTPI -Pip*- &*)* (p 2 ~ C 2 )* (c 2 - 6 3 ) O) 2 ^. 

Similarly the parts of M, N arising from this part of p are 

+ tV ^-pi • p(p^ — & 2 )* • (p 2 — c 3 )* ( — c 2 ) eo 2 2 and 0. 



(c.) Lastly, if 



P 



P\$ = — Pi 



e 



a< 



^ B i i x i - Vi - ffl} + Bjaj^ + -^a^ + -^ y^ 



6 2 



(p 3 - 6 2 ) (p 3 - c 3 ) 



f f Pp 2/1^1 do- = - -A-Trpj (c 3 - b 2 ) p (p 3 - 6 3 ) 1 (p 3 - c 3 ) 



.2\i Z 1 g*« 



,2 



J S 



p2 (^2 _ c 2) 
S 3 



- f f Pp ^ do- = + -h wV (p 3 - bj (p 3 - c 3 )* f rf» 



P 3 (p 3 - 6 2 ) 



Jf Pp y A do- = - A »Pi&V (P 2 ~ &V (P 3 - c 3 )* B^. 



Collecting the different parts, we obtain for the couples, provided there be no 
external disturbing force, 



B 



^ 



3 *iM 



L = - ■& Trpi". p (p 3 - & 3 )* (p 3 - c 3 )* (c 3 - ¥) A a>% + - 3 e 






0) 



3 0'i + 



B 



3 I „i\t 



M = + ^ ff/0l . p (p 3 - Wf (p 3 - c 3 )* c 3 . 



B 



2 „iA< 



0> S l/ 8 



M 



B, 



o) 3 ^', 



^e 



«M 



N = - tV ff/0l . p (p 3 - 6 3 )* (p 3 - c 3 )* & 3 . B x e iM 
= (p - v) B/" 



J 



> • (37), 
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where, for brevity, we have put 

v = -frir Pl .p(ffi-V)*(ffi-<*)*(c*-W) 9 ii^&w. p( P *-V*)*(p*-c*)*<* . (38). 

The terms in or arising from (b) are due to the centrifugal force of the fluid, and 
occur in consequence of the axis of rotation of the shell not accurately coinciding with 
that of the fluid, while the terms (c) are due to the effective inertia of the fluid. 



§ 8. Dynamical Equations of Motion of the Shell. 

Let A, B, C be the principal moments of inertia of the shell ; p, q, r the angular 
velocities about the principal axes. 

The position of the shell at time t + St may be found from its position at time t by 

giving it a small rotation ooSt about Oz, followed by rotations — — St, ^SL — -St 

# ^ s dt ' clt ' dt 

about the positions of the axes Ox, Oy, Oz at time t + S*. 

The direction-cosines of these latter axes referred to their position at time t are 

(1, a>8t, 0), (a>St, 1, 0), (0,0, 1). 

Hence, resolving the rotations in the directions of the axes Ox, Oy, Oz at time t. 
the component angular displacements are 

e x st, e 2 st, (a> + e 3 ) st, 

and the angular velocities about the axes Ox, Oy, Oz are 

• • • 

v x , 2 , O) + # 3 . 

Resolving these about the axes Ox x , Oy l9 Oz l we see from the scheme (35) that 

p=z0 x — co0 2 , q = 2 + o)0 l9 r = <u + ftj . . . . (39). 

Euler's equations of motion are 

Ap — (B — C) qr ■= L, 
Bq - (C - A) rp = M, 
Or — (A — B)pq = N, 

where, if the system be subject to no external disturbing force, L, M, N have the 
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values given in (37). If there be any external attracting system we must add to the 
right-hand members the couples due to this system,, estimated as though the system 
were rigid throughout. In dealing with the "free" oscillations, however, we may 
omit these terms. 

Introducing the values of p, q, r from (39) and omitting small quantities of the 
second order 

A# x - (B - C) coW l ~- c4 (A. + B - C) = L, 
m 2 + (C - A) w% + g)^ (A + B - C) = M 3 



Cft 



= N. 



Hence, on replacing L, M, N by their values (37), putting X = 9\e ik \ &c, and 
omitting the time factor, we obtain 



6\ {AX 3 + (B - C) *>*} + ih»ff 2 (A + B - C) - v \o?d\ + -*j 



B\ {BX* + (A - C) g>*} - i\a>6\ (A + B - C) - /x ] <**ff % 



e\ {ex 3 } 



(jLL — V) Bj_ 



from (32) 



Bj^ 



XW 



2o* - P 



A f 

V Q. 


















. (40) 



Hence the freedom defined by the coordinate # 3 is neutral, and a disturbance, 
which causes the coordinate 3 to vary, will not give rise to an oscillatory motion. 
From (33) the values of B 3 /t, B 3 /t are given by 



Sl ( 2p * - c '») - ^ ^ (p 3 - c'*-) - XV 3 0' s = 



X T 



Is. ( 2p 8 - 62 _ c '*) + 2 ^ ?* ( p * 



\ T 



C 



^ 



) + x 2 (c s - & s ) 0' x = 



. (41). 



J 



Eliminating ff l9 0' s , B 2 /r, B 3 /r from (40), (41) by means of a determinant, the 
periods of free oscillation are given by 



AX*+(B-C-^)a> 3 , 



^Xq>(A+B~C), 



ika> (A+B-C), BX S +(A-C-/*) a>\ 



0, 



X* (c s -6 3 ), 



XV, 



0, 



2p 



0, 



/x. 



a_^ 



^ 







X 



(P 3 ~ c ' 2 ) 



2&n 



(^-.c ,a ), 2/o a -6 a -c 



^ 



= 0. 
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§ 9. Reduction of Period Equation. 



Expanding out this determinant, we obtain 
[{AX 3 - (C - B + v) a. 3 } {BX 3 - (C - A + ft) a> 3 } - a. 3 X 3 (A + B - C) 3 ] 



X 



(2p 3 - c' 3 ) (2p 3 - W - c' 3 ) - ~ (p 2 



X 3 



c' 3 ) 



3 



+ /xc 3 X 3 {AX 3 - (0 - B + v ) w 3 } (2/ - 6 3 - c' 3 ) 
+ v (c 3 - 6 3 ) X 3 {BX 2 - (0 - A + p) w 3 } (2p 2 - c' 3 ) 
- 2w 3 X 3 (p 3 - c' 3 ) (A + B - C) {fx (c 3 - W) + w 3 } 
+ ^c 3 (c 3 - V) X* = . . 



♦ f « * 



Now since 



.2 



C 3 



4ft> ? ' 



r 3 (p 3 -e' 3 ) = 1-^- (p 3 -c' 3 ), 



X 3 



(42). 



(2p 3 - c' 3 ) (2p 3 - & 



3 



4&> 3 



c' 3 ) - ™ (P* - c*) 8 



= p 3 (p 3 - & 2 ) + (p 3 + p 3 - 6 3 ) (p 3 - c' 3 ) + (p 3 - c' 3 ) 3 (l 
= "V -WV ~ V) + (V - & 3 ) (p 3 - c 3 ) + (p 3 - c 3 ) 3 } 



1 



= ~ . {X 3 (2p 3 - c 3 ) (2p 3 - & 3 - c 3 ) - 4< y 3 p 3 (p 2 - 6 3 )} 



2p 3 — 6 ; 



3 -c' 3 



1 



.{ T 3( /3 3_&2) + (p 3 - C 3 )} 



= TO • i* (V - & 2 - c 3 ) - 4c 3 (p 3 - fc 3 )} 



XV 



and 

o„3 



2p* — 



c' 2 = ^.{rV + (p 3 - c 2 )} = ^.{X 2 (2p 3 - c 2 ) - 4c 3 p 3 } 



4« 3 ' 
X 3 



>. (43). 



Hence substituting in (42) we obtain the following cubic for X 3 : 

[( AX 2 - (C - B + v) a) 2 } {BX 2 - (C - A + p) w 2 } - ft. 2 X 3 (A + B - C) 3 ] 

X {X 3 (2p 3 - 6 3 - c 3 ) (2p 3 - c 3 ) - 4o) 3 p 3 (p 3 

+ jac 3 X 3 .{AX 2 - (C - B + v) a. 2 } {X 2 (2p 2 - ¥ - c 2 ) - 4a. 2 (p 2 - 6 2 )} 

+ v (c 2 - 6 2 ) X 2 .{BX 3 - (C - A + p) a> 2 } {X 3 (2p 3 - c 3 ) - 4ft> 3 p 3 } 

- 2aj 2 X*.(p 2 - c?) (A + B - C) {/* (c 2 - & 3 ) + vc % } 

+ p*>X 4 .(X 3 - 4ft) 2 ) c 2 .(c 2 - & 2 ) = 0. 
MDCCCXCV. A. 3 S 



& 3 )} 



J 



498 



MR. S. S. HOUGH ON THE OSCILLATIONS OF A 



Let us now change our notation, replacing p % , f? — 6 3 , p" — c 2 by a 2 , /3 2 , y 3 ; the 
period equation may then be written 

[{AX 2 + (B - C - v) «> 3 } {BX 2 + (A - C -/*) o 3 j - a> 2 X 3 (A + B - C) 3 ] 

X [X 3 (a 3 + y 2 ) (/6 3 + y 3 ) - 4a> 2 a 2 / 8 2 ] 



4. ^ ( a a _ y*) x 3 { AX 3 + (B - C 

+ „ (/3 3 - y 3 ) X 3 {BX 2 + (A - C 
- 2o> 2 Xy . (A + B - C) {/x (/3 2 ■ 

+ fxvk 4 (X 3 - 4o> 3 ) (a 3 - y 2 ) (£ 2 - 



- V ) CO 3 } {X 3 (/3 3 + f) 

-p)io*} {X 2 (a 2 + y 3 ) 
•y 3 ) + v{a 2 — y 2 )} 
y 3 ) = . . . . 



4ft> 3 y8 2 } 

4ft) 2 a 3 < 



• ( 44 )- 



If we put X 3 = or, the left-hand member of (44) becomes 



w 



6 



[(A + B - C 



,) (A + B - C - /*) - (A + B - C) 2 ] [(a 3 + y 2 ) (/S 3 + y 3 ) 
+ a>V (« 3 - y 2 ) {A + B - C - ^} {y 3 - 3y8 2 } 
+ a/V (yS 3 - y 2 ) {A + B - C - /*} {y 3 
2o)V (A + B - C) {^ (£ 2 - y 3 ) + v (a 2 



4a 2 /3 2 ] 



3 a 3 } 



y 3 )} - 3/W ; (a 2 - y 2 ) (/3 3 - y 2 ) 



oj i] ( A + B - C) > {(a 2 +/) (/3Hy 3 ) ~ 4a 3 y8 3 + (y 2 ~ a 2 ) (y 2 ~ 3yS 2 ) - 2y 2 (y 2 -£ 3 )} 

+l/ f( a 2 +y 2)( | 82 +y 2)_4 a 2 y 8 2 + (y 2 _ / 8 2 )(y 3 -3a 3 )-2y 3 (y 2 -a 2 )} 



+ jJLVOJ 



6 



+l/ (( a 2 +r 2)( | 82 +r 2)_4 a 2 y 8 2 + (y 2 _ / 8 2 )(y 3 -3a 3 )-2y 3 (y 
"(a 3 + y 2 ) (yS 3 + /) ~ 4a 3 /3 2 + (y 3 - a 2 ) (y 2 - 3^ 2 ) _ 
• + (/ - £ 3 ) (y S - 3« 2 ) - 3(y 2 - « 2 ) (y 2 - £ 2 ) 







therefore X 2 = w 2 is always one root of (44). 

Arranging (44) according to powers of X 3 , we have 



X (i 



0) 



2 X*. 



AB (a 2 + y 3 ) (/3 3 + y 3 ) + ix A (y8 3 + y 2 ) (a 2 - y 3 ) + vB (a 2 + y 2 ) (£ 2 - y 2 )" 
+ pv (a 2 - y 2 ) (yQ 3 - y 2 ) 

> 3 + y 3 )GS 2 + y 3 ) {(A + B~C) 3 -A(A-C~^)~B(B~C-v)}+4ABa 2 / 8 3 " 

- /* (« 3 - r 3 ) {(/3 3 + y 3 ) (B - c - v) - 4A^ 2 } 

- v (/3 3 - y 2 ) {(a 2 + y 2 ) (A - C - /*) - 4Ba 2 } 
+ 2y 2 (A + B - C){/t (/8* - y 3 )+ v (a 2 - y 3 )} 
+ ifiv (a 3 - y 2 ) (/3 2 - y 3 ) 



+ to*X 3 . 



~( a a 4. y») (/3 2 + y») (B - C - v) (A - C - fi) 

+ 4a 2 /3 2 {(A + B - C) 3 -A (A - C - ft) - B (B - C - v)) 

_- 4 j u/3 2 (a 2 ~ y 3 ) (B - C - v) - 4m 3 (/3 2 - y 2 ) (A - C - /*) __ 
- ft , 6 {4 a 2 / Q 2 (B-C- z ,)(A-0- j a)} 



= 0. 



Dividing out by the factor X 2 — or 1 , we are left with 
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-co 3 X 3 



X 4 . [A (/3 2 + y 3 ) + v (£ 2 - y 3 )] [B (a 2 + 7 3 ) + li (^ - y 3 )] 
"(a 2 + y 3 ) (/8 2 + y 2 ) (B - C - v) (A - C - y) 

4-4a 2 /3 2 {(A+B-C) 2 -A(A-C-/x)-B(B-C-v)-(B-C-^)(A-C~/x)} 

-4/A/3 2 . (a 2 - y 3 ) (B-C-v) - 4 w 2 (/3 2 - y 3 ) (A - C - p) 

+ o 4 . {4a 2 /3 3 (B - C - v) (A - C - /*)} = 0, 



or 



X 4 . [A (/3 2 + y 2 ) ■+ v (/8 2 - y 2 )] [B (a 2 -f y 2 ) + /x (a 2 - y 2 )] 
- co 3 X 3 ~(« 2 + y 2 ) (0 2 + y 3 ) (B - C - v) (A - C - ft) + 4a 3 £ 2 (A + v) (B + fi)" 
_+ 4/t/8V . (B - C - v) + 4m 3 y 3 (A - C - p) 

+ a>*{4 a 2 / 8 2 (B - C - v) (A -C - /*)} = (45). 

In order that the system may satisfy the criteria for " ordinary " stability the roots 
of this equation in X 3 must be real and positive. 

The period equation (45) agrees with the equation (5) (§ 1), and the solution of it, 
in the case in which the ellipsoid differs but slightly from a sphere is given in § 3. 



§ 10. Nature of the Oscillations. 

From equations (40), (41), (43) we see that the equations giving the ratios of th 
quantities 6\, 6\, B 3 /t, B s /t are 



e 



6\ {AX 3 + (B - C) w 3 } + tXo, 0\ (A + B - C) - v Lw\ + 



B 



s 



B 



B\ (B\ s + (A - C) a>*} - iXco 6\ (A + B - C) - /£ (coW, + 

\ " T 











B 



3 



T 

B 

T 



{(a 3 + y 3 )X 3 -4o) 3 a 3 } -2 W -iXy 3 - 3 -X 3 (X 3 -46i 2 )(a 2 -y 3 )^ 3 = 



- • (46). 



B, 



? {(P* + T 2 ) ^~ 4o> 3 /8 3 } + 2a)i\y* ^+ X 3 (X 3 - 4co 3 )(/3 3 - y 3 ) & - 



(a). We have seen that in every case X 3 = &> 3 is one root of the period equation ; 
when X = « the equations (46) reduce to 



0\ {A + B - C - *} + ^' 3 (A + B - C) - v 



B\ {A + B - C - p] - iff, (A + B - C) + 



/* 



B, 



TO)' 






TO)' 



0. 



o, 






B 



(y» _ 3 «3) _ 2 iy* -L + 3 (« 3 - y 3 ) 0' 2 = 0, 



B 



JSL (y3 _ 3/32) + 2 ,y 



B 



2 



TO)" 

3 S 2 



3 (/3 3 - y 3 ) 0\ = 0, 
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which are satisfied by 



p p 

TGT TftT 



Suppose ^ = <£e* M+e) where <£ is a real quantity ; we have as one solution 



l = <j>e i(bit+ *\ 2 = i<j>e 



i(<at + e) 



Changing the sign of i wherever it occurs, another solution is given by 

l = <f>e - * M+e \ d % '= - t^e " * M+e) 

and this corresponds to the value — o) of X. 

Combining these two solutions we get as the real motion corresponding to the root 
X 3 = «> 3 , 

6 X = 2<£ cos (a)£ + e), # 3 = — 2(f) sin (coi + € )* 

Now l9 0% serve to determine the position of that principal axis which in the 
steady motion coincides with the axis of rotation, relatively to axes which are 
themselves revolving with angular velocity co. 

Let us consider the angular displacements relatively to fixed axes Of, Or), 0£ 
coincident with the position occupied by the moving axes at the time t = ; they 
are 

X cos cot — 2 sin cot = 2<j) cos e, 

0] sin cot + 0% cos cot = — - 2<j> sin e, 

and these are constant quantities. Thus, the apparent oscillation which corresponds 
to the root X 3 = w 3 , consists of a small permanent displacement of the axis of rotation, 
and the system rotates as if rigid about an axis which does not accurately coincide 
with our axis Oz. It is obvious that if co be the angular velocity of rotation about 
this axis, the system and the moving axes Ox 9 Oy 9 Oz will return to their original 
positions after an interval 27t/g>, and, therefore, the system will appear to oscillate in 
a period 27r/o) relatively to these moving axes. 

It is easy to see that the fluid motion, indicated by the analysis, also consists of a 
motion of pure rotation. 

For when #, = <f>e i(oit + e) 

B x = B 2 /™ 3 = - iBJra) 2 — i<f>e u . 

Therefore, from (34), 

\Jj } = G)H<j)e ie {xz + iyz} 9 
and from (6), 
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■ u x = — ( — <jxohe ie — 2<jxohe ie ] = — <£<a2e ie 

1 

^ = ^~ { — . icf)G)he ze — 2i<jxohe le } = — cfxozie 15 

whence 

: ^ = — (j>coze i{(at + e \ v = — <j)0)zie i(oit + e) , w = <£o) (# + ^) e { ' w * + e) , 

arid, in the corresponding real solution, 

u = —> 2<f)coz cos (o)2 -{- € )> v = 2^co2; sin (cotf + e )> 

w = 2^0) [a? cos (o)t + c) — y sin (cuf + e)]. 

These are the component velocities relatively to the moving axes. The velocities 
parallel to the instantaneous positions of the moving axes are 

u — yco, v + wo*, w. 

The velocities parallel to the fixed Of, O77, 0£, are therefore 

(u — yw) cos cot — (v + ft®) sin cot = — 2<f)0)Z cos e — on?, 
(^ + xco) cos o)i + (^ — - yco) sin cot = 2(/>co£ sin € + *>£ 

w = 2c/)0) {f cos € — rj sin e}. 

Thus the fluid motion is a motion of pure rotation, the component angular velocities 
about the axes being 

— 2cf)0) sin e, — 2<j>o) cos e, a). 

The resultant of these angular velocities is an angular velocity o> about the line 
whose direction cosines are 

— 2(f) sin e, — 2<j> cos e, 1. 

The similar case for the spheroid with a free surface has been already discussed by 
Bryan (' Phil. Trans./ 1889). 

(6). Next take 

X = co [1 + ( ei +e s ) (1 + q)J = co (1 + E) say 
where 

E = i(€ 1 + e 8 )(l + (?). 

Substituting this value of X in (46) and putting a = y (1 "+ € i)> & c - 5 we obtain 
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^ 1 {A + B-C + 2AE-2C6j}+^' s {(A + B-C)(l+E)}- ? Ce 8 ^ = J 



-2 2 



^ 2 {(A + B-C)+-2BE-2Ce 1 }--i0' 1 {(A + B-C)(l+E)} + gCe 1 



B, 



TO)* 



o, 



B, 



. B 



V 



TO) 



2 



3 
Tft> 



{1 + ei ~E}+*— | (l + E)— 6»' 3 e 1 (3 + 2E) = 0, 



B a <-> ■ - .E}-iA (l + E) + ^ 1 e 2 (3 + 2E) = 0. 



TO) 



(1+e 



2 



Tft>' 



_*> 



(47). 



These equations are correct as far as first powers of e l5 € 2 , E only. Solving for the 
ratios of 0\ 9 0' s , B 2 /ro> 3 , B 3 /W 2 ? we have 



o\ 



-e' 



A + B-C + 2BE 
— oe 1; 

0, 



qGe v ^Cej, 

1 + 6 X -E, i(l + E) 
-*(1 + E), l + e 2 -E 



»(A + B-C)(1 + E), gOe 1; 

1 + ej— E, 

e 3 (3 + 2E), -i(l + E), 




i (1 + Ej 
14 « 2 — E 



B„/t&) 2 



■»(A + B-C)(1 + E), A + B 

0, 
3e 3 , 



G + 2BE~ ( /Ce I , 



3e 



i> 



0, 




■i(H-E) 
l + e„-E 



■B 3 /t«o 3 



i(A + B-0)(l + E), 

0, 
3e 2 



or 



o\ 



■% 



iff. 



A + B-C + 2BE- !Z Ce 1 , 
0, 



"kBg/T&T 



-<(1+E) 



'B 3 /to) 



2 



(A + B-~C){e 1 + e 3 -4E} (A + B-C){e 1 + e a -4E} 3(e 1 + 6 g )(A + B-C) 3(e 1 + e 2 )(A + B-C) ' 



where the denominators are correct as far as first powers of c ]9 &c, only. 



Replacing E by its value Jr (e x + e 2 ) (1 + q) we obtain 



^ 



id\ 



— iB 2 /ro) 3 



(2g +1) 2^+1 



3 



Taking 



we have 



6\ = <j>e u 



Q' % = i<f>c u 



and, in the corresponding real motion, 



6 1 = 2<£ cos (X£ •+ e), 
# 2 = — • 2<£ sin (X£ + e). 



'BJray 



3 



* • ♦ • e ITU li 



The angular displacements, relatively to the fixed axes 0£ Otj, 0£, at time £ are 
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X cos cot — 2 sin o>t = 2(j> cos [(X — a>) £ + c], 
^ sin <w£ + Q % cos <»£ = — 2cj> sin [(X — a)) £ + e]. 

Thus the motion of the principal axis consists approximately of a simple harmonic 
motion in period 2rr/(\ — o>) = 27r/Ea> 3 in virtue of which it describes a small cone 
about its mean position in a direction opposite to that in which the system is rotating. 

The position of the instantaneous axis of rotation of the shell is defined by the 
direction cosines 

or 

— 2<f> — sin (Xt + e )* — 2(p — cos (kt + € )> 1 ; 

and since X is approximately equal to co, this axis will be very nearly coincident with 
the principal axis. 

From (48) we have also 

s _ • n s , ° :a f ° 






1 + 2^ x 1 + 2q r ' 



and therefore from (34) 



o 
O 



i/f x = — 1 t<^(w a (a?2 + %^) e *" ; 

1 -f* -£j2 



to the same order of approximation we have 

o o <"> 

Mi = f ^ <M e *> V l = f^2^ ^ e "> o>i = ~ Y+: 2q ^° ^ + ^ 6 "- 

Whence, in the corresponding real solution 

u = t- ^>w2 cos (X£ + e), v = — T~r7r < ^ <uz s ' n (^ + e )> 

ty = — <£o> [ — x cos (Xi + € ) + V s * n (^ + *)]• 

The velocity components relative to the fixed axes are therefore 

6 



(u — yo)) cos cot — (v -\- #<«) sin g)£ = -~- -~- <jxoz cos (X -— a)£ + e) — a)?? 



6 



(t> + aw) cos tui + (w — yoi) sin gh5 = — 7-717- 4> MZ s ^ n (X — w£ + e) + «£, 

1 -f" ^j<7 

6 . 

w = — T"V"oi <jf>co {£cos (X — <o£ + e) — 17 sin (X — e)t -f- c)} * 

* x 
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Hence, the motion of the fluid will consist very approximately of a rotation, as if 
rigid, with angular velocity o) about the line whose direction -cosines are 

6 t m . . m s 6 



(j> sin (X — cot + c), - — — <£ cos (X —- co£ + € ) 3 !• 



1 + 2q r v ' /? 1 + 2^ 

This axis will itself describe a cone in period 27t/Eg>, and it will be so situated that 
the axis Oz lies in the plane containing the axes of revolution of the fluid and of the 
shell, and is between these two axes. The semi-vertical angles of the cones described 
by the axes of rotation of fluid and the shell will be in the ratio 3 : 2q + 1. 

The motion under discussion is that which w r ould ensue if the shell were set 
rotating about its principal axis, while the fluid possessed a rotatory motion in the 
same period about some other axis. It is clear, that as e l9 € 3 , and consequently E, 
diminish, the period of this oscillation will be prolonged ; that is to say, the motion of 
the axes of rotation will become slower. This motion will be reduced to zero when 
Cjl, e 3 vanish. In this case the internal surface of the shell is spherical, and the 
shell and fluid, of course, move independently. So far as this (apparent) oscillation 
is concerned, they will each be rotating with angular velocity co, but about 
different axes. 

From the expression for the ratio of the amplitudes, we see that when q is large, 
that is wdien the effective inertia of the fluid is large, compared with that of the shell, 
the disturbance of the shell will be considerable, compared with that of the fluid ; 
whereas if q be small, the disturbance of the shell bears to that of the fluid, a ratio 
which approximates to, but is always in excess of, 1:3. 

This oscillation has been previously examined by Hopkins ('Phil. Trans./ 1839) 
under certain special assumptions, as to the initial circumstances, and to the law of 
distribution of density in the shell. 

(c). Lastly, suppose X = cov^/q + # e i) (^2 + #e 2 ). 
The approximate form of equations (46) is now 

T> 

e \ (*2 + 1^} B - iffzB y/( Kl + qej (k^H- 2 e 3 ) + qCe. z f = ; 

&% { *) + 2 e i 1 A + i0\ A \/W+qeM^+qec,) ~ gGe, ^ = 0, 
2 ~ (1 + 2e x ) + i VT«Y+ qeJ(K,+ qe^ = 0, 
27(1 + 2e 2 ) - i V\k x + ge,) (k, + ^) ^- = 0. 

Hence approximately B 2 /r = B 3 /r = 0, and 

6\ %6' 
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therefore 

The corresponding real solution is 

0^ = 2<£ v /(ac 1 + qe x ) cos (Xt + e), 
6 % = 2(f) v/(/c 2 + ge 2 ) sin (A* + € ). 

The angular displacements about Of, O77 are therefore 

X cos cot — 2 sin cot = cf> {x/(k x + gcj) + ^/(^ 4" S^)} cos [(^ + ^) * '+ e ] 

+ <£{ \/(*i + 2 e i) — \/( K -2 + ^a)} cos [((» — X) * - e], 

2 cos o)2 + a sin o)2 = <f>{ x/(k x + 2 € i) + x/^s + 2%)} sin [(co + ^) * + € ] 

+ <£ { v/(^i + q*i) - ^/(k 2 + ge 2 )} sin [(o> - X) *. — e] 

The motion of the principal axis in space consists, therefore, of a combination of 
two simple harmonic motions, the period of each being approximately equal to 
the period of rotation of the system, and the amplitudes being in the ratio 

\/( K i + 2 e i) ^ \/( K z 4* 9 € s) : y/( K i + # € i) 4 \/( K 2 4~ # e 2) 5 * n virtue of each of these 
oscillations, the principal axis will describe a cone of revolution in the direction in 
which the system is rotating. In the event of the system being symmetrical about 
the axis of rotation k x = k% and € x = e 2 , in this case the amplitude of one of the 
oscillations reduces to zero. 



We have likewise 



1 = — 2cj> v /(ic l + qe } ) sin (Xt + e). 

ft) ft) 



+ 2(f) </(k. 2 + ^€ 2 ) cos (Xt + c). 



2 

ft) ft) 



~- cos ft>£ sin cot = — <£ {-v/^i + 7C1) + vX^s + c l € z)} sm ( w 4" ^ 4* € ) 

— <£ — Wi*i + # 6 i) - \/(^2 + 2*a)} sin (v 7 ^ - Xt - € ), 
« • » 

9 \ 

1 sin o)2 -f — cos o;2 = <j?> — { x/l^i + 2 e i) 4* \/(^2 4 <7 e 2)} cos {( w 4" M * 4 e] 

+ ^— {^/(^ + q €l ) - v/{fc 2 + £e 2 )} cos {(co - X) t - c}. 

The motion of the instantaneous axis of rotation of the shell is therefore in all 
respects similar to that of the principal axis, but the semi-vertical angles of the cones 
described are smaller in the ratio X : co. 

MPCCCXCV.— - a. 3 T 
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The direction-cosines of the instantaneous axis referred to the principal axes of the 



shell Ox ]5 Oy l9 Oz x are 



* * ; 



CO 



O) 



1, 



or 



— 2<£ 



\ 



0) 



\/( K i + 2 € i) + V'( K 2 + 9 e z) sin ( Xt + € ). 



2<j) 



Q) 



V(** + 9%) + \Z(*i + 2 € i) cos (^ + € )> i ; 



therefore, relatively to the shell, the instantaneous axis descrihes a cone in period 27r/X. 
This motion would ensue if the shell were started rotating about an axis not 
coincident with its principal axis, and it is analogous to the motion of a rigid body 
under no forces when slightly disturbed from a motion of pure rotation about a 
principal axis. 



